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Ekbata m Part I. (^Proceedings,' A, vol. 92.) 

line 23, lofite Z/T for ZP. 

equation 17, omit factor (l—x). 

equation 18, omit factor (L+1 — 2 a:). 

line 23, omit '' since |(L + 1) is greater than L." 

line 15, write clZjdt for dzjclt. 

line 21, write ^{Y-'E-r)/(K-Y)for ^(Y-^-br). 

equation 29, W7%te x = L— (L'— L)/(2/~l). 

line 22, writeY—v for v—Y. 

r2K5e~^-^^^'^ 
line 2 from bottom, torite — — -"zi^em ^^* 

equation 36, write e-^'^i for g-^icK^ 
equation 38, tvrite Po/or P (twice). 
equation 39, larite e^'^'^ for e"^^ and /or e™^^'.. 
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221, line 8, last term, lorite U for L. 

lines 4, 3 from bottom, write " seen that all the tangentials of x are 
zero." 

222, lines 6, 5 from bottom write '' dxo/dt == K(L'Pq — 1)/'Po^,^^oy nearlj 

KL/Po." 

223, line 10, turite KL/Fo for K/Po. 
225, equation 67, write (l—Jj) for (—-L). 

229, equation 80, write x/xQfor L/a?o- 

230, 3, write Part Yl for Part IV. 

Preface to Part II. 

Part I of this paper was published in the ' Proceedings,' A, vol. 92 (1916), 
having been read on ITovember 11, 1915. In June last, the Eoyal Society 
was kind enough to give a Government Grant for providing me with 
assistance in order to complete the paper, and for carrying on further 
studies upon the subject; and Miss Hilda P. Hudson, M.A., Sc.D., was 
appointed for the work from May 1, 1916. The continuation of the 
paper has accordingly been written in conjunction with her ; and I should 
like to take the opportunity to express my obligations to her for her 
valuable assistance,* especially in regard to Part III- — which is to appear 
shortly. 

I must apologise for the rather numerous small errors in Part I — due to 
the fact that the proofs were received by me when I was abroad on active 
service. 

The entire paper here presented is still limited to the theoretical side of 
the subject, as defined in the Preface to Part I. Becords of epidemics are 
now being examined in order to find how far the theoretical results which 
we have reached may be applied to them; but these studies must be reserved 
entirely for future discussion. — R. P. 

VIII. 

(i) Hypothetical Uijidemics : Deduction of Constants. — In any case of inde- 
pendent happening we can use equations 8 and 12 with V, v, 'N, r, all zero. 

Then 

dx/dt = h(l-'X), and, if Xq = 0, 

ht = — log^(l—x), 

= x-\-l'X^ + ^x^-\~ ... if X is small, 



1 

2*30 logio -— for all values of x. 

1—x 



s 2 
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If the value of x is known for any value of t, say a year, this equation 
enables us to find h, and thence x for any other period, hour, day, week, etc. 

Thus suppose that child-birth, death, or migration occurs to 20 'per mille 
of a population of all ages during one year, to what proportion do they occur 
in one day ? We have 

365/?.=: 0-02 + 0-0002 + 0-0000027+..., 
h r= 0-000055. 

To find the proportion for one hour, we divide by 24 ; and to find it for 
7 or 30 days, we multiply by these numbers respectively. 

Or suppose that 85 per cent, of plague cases, taken at any stage of the 
illness, die within 30 days, and we require the daily death-rate of such 
mixed cases. Then we have 

30/?.= -2-30 logio (1-0-85), 
h = 0-063, 

giving a constant death-rate of 6-3 per cent. jp^7' cZi^w. If the cases are all 
taken from the commencement of the illness, the curves are quite different, 
and show the highest death-rate some days after the onset — but this does 
not concern us at present. 

If the annual birth-rate is 30 per mille, and the annual death-rate is 20 
per mille, then the corresponding daily rates are 

n = 0-000083 and m = 0*000055. 

Hence if the immigration and emigration rates counterbalance each other 
w^e may take for the daily variation element v = n ■—m = 0-000028. 

In this example the increase of the population at the end of the year 
would amount to 1 per cent. — that is, P would vary from Po at the beginning 
of the year to I'Ol Po at the end of it. But if we use the value of v just 
given for equation 13, namely, P === Po^"^^ (putting ^ = 365 and 'y^ = 0-0103), 
we obtain P = 1-0103 Po at the end of the year — which is a little greater 
than the figure just given. The reason for this is, of course, that in 
equation 13 we assume that the propagation is continuous — that is, that 
during the year the progeny of the original population will themselves have 
progeny. But this is not assumed in the value of P first obtained. 

We may also calculate v from equation 13 if we know the values of Pq 
and P at the beginning and end of t. Thus the population of England and 
Wales was 10,164,256 in 1811, and 36,070,492 in 1911. Then since 
vt = logeP— logePo and' t = 365 x IDO, the daily variation-element is 
'y = 0*0000347, the yearly variation-element is 0*012666, and the ten-yearly 
variation-element 0-12666. If a population doubles itself in a period t, we 
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have vt = log^2 = 0*693, so that at the above rate the English population 

ought to double itself in 19,975 days, or nearly 55 years— which is about 
the case. 

If the daily death-rate of plague is taken at the figure just calculated, 
of 0*063, and the plague birth-rate is a normal one of 30 ^e^^ m^Z/e a year 
(so that N = n\ then V = N-M = - 0-063; so that, ii v = 0*000028, 
v—Y = 0*063.., daily. Such figures give us some concrete ideas of what 
the constants are likely to be in cases of human epidemics. [Compare 
Sections III and V (iii).] 

The reversion-element r is particularly elusive since we seldom have 
known figures regarding it to deal with. It is best to proceed as we did in 
connection with equation 20, and to assume that 90 per cent, of the affected 

individuals revert in the time t, where t = 1, 2, 3,..., 10, 'or more years. 

1 
Then for a daily element 365 r = - 2*30. Hence if 90 per cent, of the cases 

revert in one year, the daily reversion- element is r = 0*00631. Thus the 
reversion-element is likely to be large compared with 'y or N" and corre- 
spondingly important, and will be by no means negligible even when the 
considered disease is one in which the acquired immunity is supposed to be 
very lasting [compare Section III (iii)]. For a case reverts not only when it 
becomes again capable of showing a recognisable infection, but when it is 
again able to harbour the infective agents sufficiently to afford them a nidiiSr 
whence they may infect others ; and it is quite possible and indeed probable 
that this may often occur much earlier than we imagine. Thus persons 
vaccinated against smallpox may acquire the disease in a very mild and 
modified form after only a few years, and may then infect non-protected 
persons in full force ; and there is no proof that a mild first attack of 
measles or other diseases may not give quite as short-lived an immunity — 
that is, that those who have had the disease in infancy may not compara- 
tively early acquire it again in an indistinguishable form and then spread it. 
In such cases, r may be considerably larger than we might otherwise expect. 

Eegarding immigration and emigration, it will suffice to note that in 1911 
there were 350,429 immigrants into the United Kingdom and 454,527 
emigrants from it. 

Further remarks on the constants will be found in Section X. 

(ii) We may now proceed to consider numerical examples of proportional 
happenings. For this purpose it is convenient to adopt the suggestion made 
towards the end of Section YII (vi), and to take c as a function of a,n 
independent parameter 7 which is such that when 7 is positive KL is 
positive and x always increases with the time, and when 7 is negative KL is 
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negative and x always diminishes with the time (Section IX). We shall 
also write D for 7;--V and E for N" + r. The principal constants of x and/ 
may then be written 

c = D + (7 + l)R, K = (7 + l)Ps KL = 7E, 

L- ^ I- c "^ . „(7+iy 

7 + 1 (7 + If 47 

The reader is reminded that KL is the coefficient of / in x, the proportion 
of affected individuals (equation 50), and that the .-r-curve reaches its centre 
of symmetry when x = |L, and its maximum 'value L when t is very large 
[Section VII (iii)]. The period at which x reaches its centre of symmetry 

1 T 

and dx/dt reaches its maximum is t = r — ~— log^ -. If ^ == 1/Po and 

JAsJu Xo 

LPo is large compared with unity, then, as suggested in connection with 
equation 58, we may now write 

7ET = logePo + loge7~-loge(7+l), 

and at twice this period the maximum number of affected individuals less 
one is reached. 

Similarly, when t is very large, the ultimate value of the proportion of 
new cases is /. But if L > |-, /rises to a previous maximum value ^c when 
^ = •!•, and thereafter falls to its u.ltimate value I [Section VII (v)]. The 
ratio between the maximum and the ultimate value is \ (equation 67). The 
period at which the maximum value is reached is 

"^ ~KL ^^^^2L^ 
(equation Q6) ; and when Xq = 1/Po and LPo is large, this may be written 

7Et' = logePo + loge7-~loge(7 — 1). 

But this maximum of /(previous to its ultimate value) occurs only if 7> 1. 

It is evident that when c is expressed as a function of 7 in this manner, 
K, KL, L, T, T, and \ are independent of v and V, and that I contains them 
only as parameters of c. 

On the other hand, r and r are functions of the original population Po and 
increase as it increases — because, obviously, x and / require a longer period 
to reach their corresponding values if the population is large. It is 
convenient therefore to adopt a definite figure for Po in the following Table ; 
and we put Po = 10,000, which is, say, the population of an ordinary town. 
We may also assume that unity is negligible compared with LPo unless 

7<0-l 

If Po is not 10,000, but some multiple or fraction of it, say 10,000^, all 
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we have to do to find t and t' is to add to or subtract from the figures 
given in the corresponding columns of the Table the quantity log^p/yB,. 
Thus, if Po = 1000, we subtract 2'30/y'R. As 7 increases above 100, 
T and T approximate more and more closely to the value logePo/7E, 
T being always less than this and r greater than it. 

The actual curves depend largely on the values of E. If we suppose the 
case of a human disease, of which the affected birth-rate E" equals the 
natural birth-rate n, which is taken at 0*000083 daily, corresponding^ 
according to the previous sub-section, with an annual birth-rate of 30 pe?"- 
mille — then E equals the former figure if r is so small as to be negligible. 
But if r is a daily rate corresponding to an annual reversion of 90 per cent.^ 
it will amount to as much as 0*0063, so that E = 0*0064. We will suppose 
that these two values of E are a minimum and a maximum. 

Now if we suppose also that c = 1, that is, that each affected individual 
daily infects or reinfects one other individual, then we have 7 = (1— •D)/E— 1. 
That is, if D = 0, 7 = 12,000 if E is a minimum, and 7 ~ 155*5 (or roughly 
150) if E is a maximum. Thus, if the infection rate is near unity, the value 
of 7 will reach from three to five figures, and the epidemic will resemble 
the last examples given in the Table. 

In Example 1, if the population remains the same during the whole of 
this long period (which is, of course, never likely to occur), only about 
ten living persons out of the 10,000 will be found unaffected at the end of it ; 
and the ultimate and maximum daily number of new cases /Pq reaches only 
0*00083, that is, only one new case in 3*3 years. 

If, however, E has its maximum value, there is one new case in about 
sixteen days. 

In Example 4, half the population is ultimately affected, but still / has no 
maximum before its ultimate value. When, however, in Example 5 
7 increases above unity, this maximum begins to appear, and the ratio X 
between the maximum and ultimate values of / constantly increases with 7. 
At the same time, r and t diminish and approach each other. 

If c is nearly unity, r = 2*30 logio Po = 9*2 days, and / reaches its 
maximum in a short time, whatever the value of E. More generally, if 
D = 0, c =: (7-1- 1)E, and if 7 is large (say over 60), 7-1-1 = 7 roughly, and 
T = 1/c . logePo- In other words, if 7 is considerable, not only does L = 1 
nearly, but the maximum {\c) of / and the mode, or time {t) when that- 
maximum is attained, depend roughly on c only, the population being fixed. 

If L>-|- the /-curve slopes downward from its summit more slowly than 
it rose towards that summit. Put x ~ ^-\-^, then / = c (I— -P), and has the 
same value for equal and opposite values of f. 
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But 



df 

dt 



= 2c 



f f-S = 7^-^<l-^)(i + ^)' 



and is less when ^ is positive than when ^ has the equal and opposite value. 
Hence, at two points of the /-curve whose ordinates are equal, the descending 
side is less steep than the ascending side. 

Since / =— — "■ + c(l — L)^, and the first term on the right gives a 
y\. dt 

perfectly symmetrical bell-shaped curve, we may call c(l — L)<23 the excess of 
the /-curve. When D = and 7 is considerable, the excess varies roughly 
from -JE at the centre of symmetry to K, when t is very large — that is^ 
roughly, from \l to I, Both of these become small when r, the reversion- 
rate, is small. 

The important result, therefore, follows that, if also c is near unity, and 
therefore 7 is large — that is, if the infection-rate is high while the reversion- 
rate is low — then the /-curve becomes a nearly symmetrical bell-shaped 
curve. This, according to Dr. Brownlee (Section I), is just the kind of 
curve presented by epidemics of certain zymotic diseases. The pmna facie 
inference is therefore that such epidemics may be wholly or chiefly mere 
cases of proportional happening, as defined in Section YII (i). It is even 
the case that, in such diseases, the fall of the curve is often more slow than 
its rise, as happens with the /-curve. 

When 7 is large and L = 1 nearly, the two changes of curvature of the 
curve occur, according to equation 65, when Qx = 3^ ^yS, that is, when 
(equation 53) 

j'Rt = logePo + loge(l/a?— 1), 

= log,Po~l-32 and log.Po + l'OS, 

the latter being nearer to the summit than the former. 

(iii) Several methods may be employed for calculating the curves of 
^,/, Z and F in detail. The simplest is to divide the period during which x 
varies from Xq to JL into a number (say 10) of equal parts and then to- 
calculate the ordinate of the curve at each section. In equation 58 we used 
T to express this interval when Xq = 1/Po, but it may be employed more- 
generally to express the abscissa of the centre of symmetry, when xo has any 
value, so that KLt = log^(L/^o — 1). We now put t = tT/10, give to T the 
successive values 1, 2, 3, ,_, 10, and calculate the corresponding values of x„ 
namely xi, x^, xz, ..., and so on. Thus we have from equations 50 

^T = L{H-7(L/^o-l)^^-^}-^ 
The root can be easily evaluated by means of logarithms, and the value of 
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X then obtained from a table of reciprocals {e.g., as in Barlow's Tables). 
Owing to the symmetry of the 5ti-curve, its values when T lies between 10 
and 20 can be easily found — for example .^13 = L— 0:^7. When we have 
ascertained the ordinates of x, we can quickly calculate those of / from 
equation 61, and, unless the case-mortality is high, we can generally assume 
for short and sharp epidemics that the original population remains constant 
and that Z = ^Po and E =/Po. It will be observed that, by this method, 
when T has been first calculated we need subsequently deal only with the 
values of L and of L/^o — 1. Thus suppose that L = 1 nearly, that P = 10,000, 
and that ^0 = 1/Po, then we have approximately, 
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for example /x =/i3 = cxriXxz- 

Another useful method for calculating x and / is to transfer the origin to 
the centre of symmetry (equation 56). But if the ordinates are required for 
successive natural units of time, such as days or weeks, in order to compare 
them with statistics of epidemics in which such units have been used, then 
we must of course give to t the successive values 1, 2, 3, ..., and calculate the 
ordinates directly from equations 50. Here the symmetry of x will not be so 
useful, because the centre of symmetry is not likely to coincide with any one 
of the periods of time for which the ordinates have been found. It will, 
therefore, be generally easier to plot the curve by the method first indicated, 
to measure off the natural units of time along the time-axis, and then to 
estimate the corresponding ordinates geometrically. 

(i) Hypometric Happening : KL Negative. — We have hitherto always 
supposed that the happening-element is constant during the whole of the 
considered period, but it is now necessary to examine the changes which will 
be produced in the functions by certain changes in this element. We here 
study only the cases of dependent happenings already discussed, and these 
only for hypometric happening. It is no longer very useful to take c as a 
function of 7. 

In Section YII (vi), while considering the constants K, L, c, we saw from 
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equations 50 that if KL is positive, x always increases and is therefore always 
greater than %, its value when ^ = 0. This implies that c > D + E. When 
this inequality becomes an equality, KL vanishes and dx^jdi = 0, that is, x 
remains near its original value, when t is not very large, so that the 
happening may be said to be isometric [see (v) below]. When, however, 
.c<^D + E, KL becomes negative; and it is seen from equations 50 that, as 
g-KL« \;^Q-^ always increases with the time, x always diminishes, and must 
therefore be always less than Xq, 

We call such cases hypometric hafp'pening. That is, we now suppose that, 
after having been greater than D -f E, c becomes less than it, or even zero. 
What, then, will be the nature of the curves ? 

The solutions of equations 44, 45, 47, 51, 52 remain valid for all real 
values of c, and we have therefore only to interpret results already known. 



KL 




It is first necessary to examine K, KL, and L more particularly as 
functions of c, but we shall do so only in the case when D is positive 
(injurious happenings') or zero. Then D, E are positive constants. We have 

K = c— D, KL = c—D— E, 

represented by parallel straight lines meeting the axis of c at an angle of 45° 
•on the positive side of the origin, 

L = l--.E/(c~-D), or (L-l)(c-D) = -E, 

and the graph of L is a rectangular hyperbola, with asymptotes L = 1, 
.c = D, dotted lines in the figure, parallel to the axes of co-ordinates. 

The last portion of this L curve, for c >>D + E, gives cases of hypermetric 
happening already studied ; but we have to consider the range = c < D -f-E : 
.since L becomes infinite when c = D, we have to consider the two portions : 
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= c <^ D, for which L is positive and greater than 1 ; and D <; c <; D + R, 
for which L is negative. 

(ii) First let e = ; that is, suppose that the happening suddenly ceases 
altogether, while all the other elements remain the same. It is advisable to 
use marked letters for the quantities after this event. We take it then that 
when the original happening c ceases, f ~ 0, and that the values which 
X, Z, P, had reached up to this time are now expressed by Xo\ Zo^ Po^ while 
x\ Z\ P' denote these functions as f increases. Similarly, K^ and L' are 
now the values of K and L without the element c. Then, since K^ = —D, 
equations 51 and 62 become 

r = Po'€^^'|l-^(l--e-"(»+i^)OT, (81) 

Z' =: Zo' 6(^~^)^'. (82). 

As V— E ~ — -M + I— -E — r, Z' decreases from its original value Zq' or Z, 
and ultimately approaches zero— unless the immigration is greater than the 
sum of the other elements (as may happen in a war area)— and its diminution 
is more rapid when M, E, r are large. With regard to equation 81, we 
observe that Po'e""^' is the natural increase of the population, living when the 
happening ceased, after that event ; and the equation shows that P^ equals 
this quantity less a function of f which is zero when f = 0, which 
increases with f, and which finally slowly approaches the limit ^oVl-'' • Po'e^^' 
when f is large. That is to say, the population which remained when the 
happening ceased continues to suffer from its effects for some time after- 
wards, and then finally reaches the value 

Here Xq^/U is always less than unity, for xq' is the value of x at the moment 
when the happening ceased, and this could not have been greater than L,. 
which is less than unity [Section VII (vi)], while, as shown in the previous 
sub-section, V is greater than unity when c = 0. Hence P' is always less 
than Po'e^'^', and the difference gives the loss of population due to the .cases 
which continue after the happening ceased. 

The fraction V/xq^ occurs in the denominator of x' in equation 60.^ As 
it is greater than u.nity, the coefficient of e~^^'^'*' is positive, so that x' 
always diminishes as f increases — as could have been already inferred from 
equations 81 and 82. 

(iii) The next case is when c' is some quantity between zero and D (which 
is generally very small). Eor the study of this case, the first of equations 50 
may be written 

X' = xo- |^ + (l-'^)6-K'i^'''y', (83) 
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which diminishes indefinitely as t' increases. If c' = D and L' is con- 
sequently infinite, this becomes x^^'^'^', which vanishes. 

The values of P' and Z' are easily obtained from equations 51 and b2^ 
if we remember that the indices D/K' and c' jK! are now negative. When 
o' = D, the expressions for P' and Z' become indeterminate, but. can be 
easily evaluated, and we have 

p/ ^ p^ ^vt'-q^ ^here ^ = Xq' D (1 - e-K*')/K. 

(iv) When c' lies between D and D + E, K^ is positive and L' is negative, 
and varies from -co to zero. Then it will be seen from equation 83 
that x' is still positive, and, as in the previous cases, diminishes indefinitely 
as t' increases. For it may be written 

where K'L' is negative. For P^ we have 

F = Po' e^*' |l-T> (l-e^'^'oV^^^'^ (85) 

which is also always positive since JJ is negative, and increases without 
limit when f is large. P' is less than Po'e^*', because the effects of the 
original happening continue to be felt by the population after the happening 
itself has been reduced from c to g\ 

(v) When c' = D + E exactly, JJ and KL' = 0, and, when f is small, 
X — xq^ nearly, and the happening may be said to be isometric, . The values 
of x^ and P' become indeterminate, but can be ascertained by finding the 
limits when L' = of the expressions in equations 84 and 85, so that 



x' = xo' ] 6-^'^'''^'^)6-^'^'^'--l) \ ^ (84) 



F ^^-pje^^'il-i-xJKVy^lK 

Or the same equations can be obtained by integrating equation 47 after 
putting L = in it, and then integrating both sides of equation 44. In 
this case also, then, x^ tends to decrease as f increases, and finally approaches 
zero, when t^ is very large. 

(vi) If K'L' is positive, however small it may be, x' always increases, as 
already seen in Section VII (iii). 

It must be remembered that, though x', the proportion of affected 
individuals, diminishes, this does not mean that new cases cease to occur. 
On the contrary, unless c is absolute zero, / is the same function of ^' as 
before. And ra^' expresses, as before, the reversions of the old cases. 

When D = 0, the phase described in (iii) above is suppressed, and we 
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begin, when c = 0, at the point when K^ = 0, W\J — —E, and \I ~ —00^, 
and the functions will be easily understood. 

(i) Parameter Analysis.- — We have now considered x, /, and P as functions 
of the time, but it remains to examine how they vary if we give different values 
to the parameters li, r, v, Y, n, m, i, e, N, M, I, E. The symbol I has been 
used to denote the limit of / when t is very large. This is an important 
quantity in disease happenings because it gives the number of new cases 
which continue to occur when x has reached its limit, that is, after epidemic 
manifestations have ceased. It is, by hypothesis, these final new cases which 
keep the infection alive permanently in the population ; we may call I the 
endemic ratio. First consider ^ as a function of r 

^~ = ~ (2L™1) = A (C-D-2R). 

or K K^ 

Thus I is positive when r = (if KL>0) and increases with r until this 
reaches a value which makes L = | and Z = Jc [Section VII (vi) ]. After 
this I diminishes as r increases. But if c— D > 2E then r never reaches this 
value and I always increases. Thus I considered as a function of r has forms 
similar to those of / considered as a function of t. 

The behaviour of Z as a function of c is of the same nature. A similar 
procedure applies to / and P. 

The effect upon the total population P due to variations in the parameters 
is important. The fundamental equation is 44 : 

d^/dt = 'yP-D^P. 

As P and ix? are always positive, the sign of dV/dt depends on v and D. If 
D = 0, the happening is equivariant and the natural change of population 
due to V is not affected. If D is negative, the happening is beneficial and the 
natural change of population is augmented by it. If D is positive, the 
happening is injurious and the natural change is reduced. If v—Dx is 
negative, especially when ^ ~ L, the total population will diminish indefinitely. 
It is found that 9P/3r is always positive ; that is, in infectious diseases, a 
quicker loss of affectedness, including loss of immunity, is beneficial. This 
may seem surprising until we reflect that, as 'y > Y, the longer the affectedness 
continues in the individual the greater will be the loss of life due to it, since 
we generally suppose in such cases that M>m and that these elements act 
during, not one, but every time-unit lived by the affected part of the 
population. Of course, the greater r is the more quickly on the average will 
affectedness, including immunity, be lost. 
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Next, dV/dc is zero or negative, and the higher the infection rate the more 
injurious is the happening. 

9P/3V is zero when ^ = 0. It generally increases at first with t and 
afterwards decreases; it is negative when ^ = L provided 1 — c(l-— L) /K is 
negative. In other words, in these cases an increase of Y is beneficial at the 
beginning of the outbreak, but injurious later on. Now V = IsT — M + I — E 
and diminishes if the case-mortality M increases ; in fact 3P/3M = — 9P/3V. 
Thus an increase in the case-mortality is injurious at first, but may be 
ultimately leneficial, because an affected individual who dies ceases to be 
infective. 
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